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ASYMPTOTIC BEHAVIORS OF FUNDAMENTAL SOLUTION 
AND ITS DERIVATIVES RELATED TO SPACE-TIME 
FRACTIONAL DIFFERENTIAL EQUATIONS 


KYEONG-HUN KIM AND SUNGBIN LIM 


Abstract. Let p(t, x) be the fundamental solution to the problem 
d^u = —(—A)^u, a G (0,2), j8 G (0,oo). 

In this paper we provide the asymptotic behaviors and sharp upper bounds of 
p(t, x) and its space and time fractional derivatives 

Dlllp{t,x), \/n&Z+, 7G [0,/3], [0,oo), 

where D" is a partial derivative of order n with respect to x, (—Aj;)^ is a 
fractional Laplace operator and and if are Riemann-Liouville fractional 
derivative and integral respectively. 


1. Introduction 

Let a G (0, 2), /3 G (0, oo) and p{t, x) be the fundamental solution to the space- 
time fractional equation 

dfu = A^u, (t, x) G (0, oo) X (1.1) 

with u(0,x) = uo(x) (and dtu(0,x) = 0 if a > 1). Here 9“ denotes Caputo frac¬ 
tional derivative and := — (—A)^ is the fractional Laplacian. The fractional 
time derivative of order a G (0,1) can be used to model the anomalous diffusion 
exhibiting subdiffusive behavior, due to particle sticking and trapping phenom¬ 
ena (see [251 HE])? and the fractional spatial derivative describes long range jumps 
of particles. The fractional wave equation cl“it = Au with a G (1,2) governs 
the propagation of mechanical diffusive waves in viscoelastic media (see |2S1133] )■ 
Equation (HI) has been an important topic in the mathematical physics related to 
non-Markovian diffusion processes with a memory [261127], in the probability the¬ 
ory related to jump processes [ 3111123 ] and in the theory of differential equations 
[ZllllITllSIllSl. 

The aim of ths paper is to present rigorous and self-contained exposition of 
fundamental solution p(t,x). More precisely, we provide asymptotic behaviors and 
upper bound of 

Dl{-A^)^Dflfp{t,x), V n G Z+, 7 G [0,/3], cr,(5 G [0,oo), (1.2) 

where // and Df denotes the Riemann-Liouville fractional integral and derivative 
respectively. 
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There are certainly considerable works dealing with explicit formula for the fun¬ 
damental solutions and their asymptotic behaviors (see e.g. m [g [ini HU na na 
MM)- However, only few of them cover the derivative estimates of the funda¬ 
mental solutions. In [i US HU HO], upper bounds of (11.21) were obtained for /3 = 1, 
a = 1 — a, and 7 = 0. Also, in m asymptotic behavior for the case a = 1,P G ( 0 , 1 ) 
and cr = 0 was obtained. Note that it is assumed that either a = 1 or /I = 1 in 
[S [ig [H [m eq], and moreover spatial fractional derivative and time frac¬ 
tional derivative D^p are not obtained in [gHgHZiEg and m respectively. Our 
result substantially improves these results because we only assume a G ( 0 , 2 ) and 
(3 G (0,oo) and we provide two sides estimates of both space and time fractional 
derivatives of arbitrary order. 

Our approach relies also on the properties of some special functions including 
the Fox H functions. However, unlike most of previous works, we do not use the 
method of Mellin, Laplace and inverse Fourier transforms which, due to the non¬ 
exponential decay at infinity of the Fox H functions, require restrictions on the 
space dimension d and other parameters /3, 5, 7 , and a (see Remark 15.21) . 

Below we give two main applications of our results. First, consider the non- 
homogeneous fractional evolution equation 

dfU = A^u + f, (t,x) G {0,oo) xM.‘^ (1.3) 


with u{0,x) = 0 and additionally dtu{0,x) = 0 if a > 1. One can show (see e.g. 
[gHllEol) that the solution to the problem is given by 

ct 

/ ^ /. ^ I./ ^ -ig 


Gf{t,x)=[ [ Qa{t-s,x-y)f{s,y)dyds 
Jo Jr'^ 


where 


Dl °‘p{t,x) : a G (0,1) 


ir"pit,x) : 06 ( 1 , 2 ). 


Qa (tj ^) ■ — 

It turns out (see |181119] H that to obtain the Lq(Lp)-estimate 

ll^^^•^llLg((0,T),Lp(R<^)) — -^ll/llL<,((0,T),Lp(R<i))j P, (7 > 1 
it is sufficient to show that for a,b > 0 


0 J\x\>b 


\A^Qa{t,x)\dxdt < N 


7,2/3’ 


(1.4) 


N 1'°° f a 

dtA^Qa(t,x)\dxdy < —, / / \'\/xA^Qa{t,x)\dxdt < Na~'^. 

One can use our estimates to prove the above three inequalities, and therefore (El 
can be obtained as a corollary. Our second application is the Lp-theory of the 
stochastic partial differential equations of the type 

d^u = Au + f g{s,x)dWs, 

Jo 

where cr < | and Wt is a Wiener process defined on a probability space {iJ,dP). 
One can show (see 0) that the solution to this problem is given by the formula 
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Here Pa- {t, x) is defined as 



:cr<0 
x) : a > 0. 


As has been shown for the case a = 1 (see [m Hi [211135]), sharp estimates of 
D^A'^Pa can be used to obtain Lp-estimate 


IIA^nlL p(nx(0.T)xR‘^) — NMl p(nx(0,T)xR‘^)- 


(1.5) 


The detail of (11.51) will be given for 7 < (2 A A^) in a subsequent paper. 

The rest of the article is organized as follows. In Section |2] we state our main 
results, Theorems \TT[\m and HD In Section 13] we present the definition of the 
Fox H functions and their several properties. For the convenience of the reader, 
we repeat the relevant material and demonstration in M and HZ], thus making 
our exposition self-contained. Section jTj contains asymptotic behaviors at zero and 
infinity of the Fox H function. In Section [5] we present explicit representation of 
fundamental solutions and their fractional and classical derivatives. Finally, in 
Section | 6 | we prove our main results. 

We finish the introduction with some notion used in this article. We write f < g 
for |a;| < S (resp. |a;| > S) if there exists a positive constant C independent of x 
such that f{x) < Cg{x) for \x\ < 5 (resp. |a;| > S), and f ^ g ioi \x\ < <5 (resp. 
\x\ > 6) ii f < g < f for |a;| < 6 (resp. |a;| > 6). We say / ^ g as |a;| —>■ 0 (resp. 
|a;| —>■ 00 ) if there exists e G (0,1) such that / ~ g for \x\ < £ (resp. \x\ > e“^). We 
write f{x) = 0{g{\x\)) as |a;| —>■ 0 (resp. |a;| —>■ 00 ) if there exists 5 > 0 such that 
\f{x)\ < | 5 (|a;|)| for |a;| < 6 (resp. |a;| > 6). We use to denote a definition. As 
usual stands for the Euclidean space of points x = {x^, ..., x"^), Rq := R'^ \ {0}, 
and Z+ := {0,1,2,---}. For multi-indices a = {ai,... ,ad) G Z'^, n G N and 
functions u(x) we set 



[aj is the biggest integer which is less than or equal to a. By P we denote the 
d-dimensional Fourier transform, that is. 


:= [ e-'(“’«)/(x)dx. 


For a complex number z, 5R[z] and 3[z] are the real part and imaginary part of z 
respectively. 


2. Main results 


We first introduce some definitions related to the fractional calculus. Let /3 > 0. 


For a function u G Li(R'^), we write A^u = / if there exists a function / G Li(R'^) 
such that 




For u G Li((0, T)), the Riemann-Liouville fractional integral of the order a G (0, 00 ) 
is defined as 




One can easily check 


( 2 . 1 ) 
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Let n S N and n — 1 < a < n. The Riemann-Liouville fractional derivative Df 
and the Caputo fractional derivative are defined as 



By definition (12.21) for any a > 0 and u € Li((0,T)), 

= u. (2.3) 


Using (I2.1I) - (I2.3I) . one can check 


d?u = D? 



n — 1 < a < n. 


Thus Dfu = dfu if u{0) = = • • • = = 0. For more information on 

the fractional derivatives, we refer the reader to |291132) . 

For (T G R we define Riemann-Liouville fractional operator D“ as 


: a > 

/]"' : a < 0 


Then by (12.11) and (12.21) . for any a, /3 > 0 

zif/f = ror^- 


The Mittag-Leffler function Ea,p{z) is defined as 

OO ^ 

EaA^) ■='^ rf L-l /^V a>/3eC,5RH>0 (2.4) 

^r(afc + /3) 

for z G C and we write Ea{z) = Ea^i{z) for short. Using the equality 

VaGM,a>0 

1 (a -I- 1 — cr) 

one can check that for t > 0 


ID^£;„(-At“) = t-^Ea,i-a{-Xt°‘), cr G M, A > 0, 
and that for any constant A, 

At) := EA-Xn 

satisfies (p(0) = 1 (also </?^(0) = 0 if a > 1) and 

d^ip = —X(p, t > 0. 

Let a G (0, 2) and f3 G (0, oo). By taking the Fourier transform to the equation 
= A^u, t > 0, u{0,x) = uo{x), (and m'(0, x) = 0 if a>l) 

one can formally get 

E {u{t,-)} = {uo} . 

Therefore, to obtain the fundamental solution, it is needed to find an integrable 
function p(t, x) G Li{R‘^) satisfying 

E{p{t,-)} = EA-\^An- 


(2.5) 
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Denote 




In the following theorems we give the asymptotic behaviors of x) as 

M —>■ 0 and M —>■ oo. We also provide upper bounds when M < 1 and M > 1. 
Firstly, we consider the case M —^ oo. 

Theorem 2.1. Let a S (0,2), /3 G (0,oo), 7 S [0,oo), tr G R and n G N. There 
exists a function pit, ■) G Li(R'^) satisfying (j2.5l) . Furthermore, the following as¬ 
ymptotic behaviors hold as M —>■ 00 ; 

(i) If P G N, then for some constant c > 0 depending only on d, n, a, /3, a 



and 



( 2 . 6 ) 


(ii) If a = 1, j3 ^ N, and a = 0, 



(2.7) 


and 



(Hi) IfjG (0,/3) \ N, 




and 



(iv) If (3 ^ N and 7 G [0, /?) fl Z+, 

|A^D^p(t,a;)| - |a:|-d-27-2/3^-a+a 


( 2 . 8 ) 


and 


\D:Amfpit,x)\<\x\ 


— d— 27 —2/3—n^—cr+CK 


(v) If ^ = (3 and d>2, 



(2.9) 


and 



( 2 . 10 ) 


(vi) If 'y = (3, CT + a G N, and d = 1, then (12.91) and ( 12 . 101 ) hold. 
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Remark 2.2. (i) Note that D'^'Dfp{t,x) has exponential decay as M —>■ oo only 
when /3 is a positive integer. 

(ii) Let X 0. Then by (12.6|) and (12.81) with 7 = 0, x) —?> 0 as t —>• 0 if 

either /3 is a positive integer or cr < a. 

(iii) If 7 = /3 and d = 1, then we additionally assumed cr + a G N. Without this 
extra condition we had a trouble in using Fubini’s theorem in our proof. 

(iv) Note that we have only upper bounds of x) unless n = 0. This 

is because, for instance, Dip{t,x) is of type x^g{t,x) (see (15.11) 1 and becomes zero 
if X® = 0. Hence we can not have positive lower bound of Dip(t, x) for such x. 


Secondly, we consider the case M —^ 0. 


Theorem 2.3. Let a,/3,j,a be given as in Theorem \2.1\ and n £ N. Then the 
following asymptotic behaviors hold as M —>■ 0.' 

(i) If"f£ [0, /3) and ct + a ^ N, 

:-{d+2~,) 


w 

|A^]D)^p(t,x)|-- (1 + |ln|xp^t““|) 

, |_d_2T.+2/3^-a-a 


7</3-f 

7 = ^-1 (2.11) 
7>/3-f 


and 


:t{d + 2y+2) 
20 


\D2Amfp{t,x)\ < <( |x|2-"r‘^-“(l +|ln|x|2/5i-“|) 

l^-j—d-27+2/3—n^-CT-a 

(ii) If ^ £ [0, 0) and cr + a G N, 

\Amfp{t,x)\ - <( |a;|-d- 27 + 4 / 3 ^-a- 2 a (l + | In |) 


\ — d— 


27 + 4 / 9 ^—cr— 2 Q: 


y = /3-|-l (2.12) 


7 < 2^ - i 
^ = 20-1 (2.13) 

7 > 2/3 - i 


and 


y{d+2-/ + 2} 
20 


|^| 2 -n^-cr-- 

|i:>!;(AT'D^p(t,x)| < { |xp-"t-'^-2“(l + |ln|x|2/3t-“|) 
|^|—d—27+4 ^—(T—2a 

(Hi) If a = 1 and a = 0, 


7 < 2/3 - I - 1 
7 = 2/3-f-1 (2.14) 

7 >2/3-1-1. 


d + 2y 


d + 27+2 


|A>(t,x)| ^ t-^, |17^A>(t,x)| < |x|^-”t-^^. 

(iv) If "f = 0 and d>2, 


and 


't-^ 




d 

2 

<0 

|x|- 

-d+2/3^- 

cr 2a -I- |ln |xp^t “ 1) 


d 

2 

= 0 


-d~\-20^— 

(7—2a 


d 

2 

>0 

[lx 

2-n^-<T 

o(d+2) 

“ 20 

d 

2 

+ 

l</3 

< |x 

2-n^-a 

-2a(i _|_ |^|2/3^-ap 

d 

2 

+ 

1 = 0 

[|x 

— d+2/3 — 

n^—(7—2a 

d 

2 

+ 

l>/3 


(2.15) 


(2.16) 














ASYMPTOTIC BEHAVIORS OF FUNDAMENTAL SOLUTIONS 


7 


(v) If 'j = P, tr + a G N, and d — 1, then (j2.15ll and (12.1611 hold. 

Next we give the upper estimates when M > 1 and M < 1. 

Theorem 2.4. (i) Assertions (i)-(iv) of Theorem \2.1\ also hold for M > 1 i/ 
is replaced by “< 

(ii) Assertions (i)-(v) of Theorem \2.S\ also hold for M. < 1 if is replaced by 

n 


The proofs of Theorems 12.1112.31 and 12.41 are given in Section [51 
Remark 2.5. Let a = 1 and fi £ (0,oo), and cr = 0. Then by Theorem 12.41 


\D:A^p{t,x)\< 


: 7 e Z+ 


I —d—2'y — n 


: 7 ^ Z+ 

holds for \x\^ > t. Also, by Theorem 12.41 for \x\'^ < t, 

^4-2-7 o d+2'Y-f2 

\A'^pit,x)\ < |Zl”AXt,x)| < |xp-"t 2^. 

These estimates cover the results of |151 Lemma 3.1, 3.3] and |131 Corollary 1]. 

Remark 2.6. Theorem 12.41 also implies the results of [HI Proposition 5.1, 5.2]. Let 
/3 = 1 and take a jnj > 1. For jxp > t“. 

\D'f^W^p{t,x)\ < |xr^-"t-'^exp{-ct-5^|x|^}. 


Q(d + Tl) 

<t 2 °’eXP'{-Ct 2-a.|a;|2-a 


'{■ 


}■ 


Also, (cf. (lOTTl . ^rm . and (ITTdll L 


t 2 

\p{t,x)\ < <( t““(l + |ln|xpt““|) 

,\-d+2t-a 


d > 3 
d = 2 
d=l 


and 


\Dlp{t,x)\ < |x 


— d-\-2—n±—OL 


, |XAp(t,x)| <|x|-‘^+2-"X“ 

I — rf+4 —n^—ct— 1 


\D'fIi]-<^p{t,x)\ < { |x|2-"t-i(l + |ln|x|2t-“|) 
|x|i-"t-i 


d> 3 
d = 2 
d=l 


hold for IxP < t°‘. 


3. The Fox H function 

In this section, we introduce the definition and some properties of the Fox H 
function. We refer to [14] for further information. 
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3.1. Definition. Let r( 2 ;) denote the gamma function which can be defined (see 
[2l Section 1.1]) for z G C \ {0, —1, —2,...} as 


ra->oo z{z + 1) ■■■ {z + n) 

Note that r(z) is a meromorphic function with simple poles at the nonpositive 
integers. From the definition, for z G C \ {0, —1, —2,...}, 


zr(z) = r(z + 1 ), 

and it holds that 

m — 1 j 

r(i - z)r(z) = TT r(z + -) = ( 2 ^)'^m^-™"r(mz). 

Sin 7 r 2 : m 


One can easily check that for k G 

Res 2 ^_fc[r( 2 )] = lim {z-\-k)T{z) 


(-1)^ 

k\ 


(3.1) 


(3.2) 


where ReSj^j;^ [/(z)] denotes the residue of f{z) at z = zq. From the Stirling’s 
approximation 

F(z) - |z|^oo, 

it follows that 


|F(a + i 6 )| - |a| ^ oo (3.3) 

and 

|r(a + i6)| - |6| ^ oo. (3.4) 

Let m, n, v, fi be fixed integers satisfying 0<m<fj,, 0<n<v. Assume 
that real parameters Ci,..., Cj^, 5i,..., 0 ^ and positive real parameters 71 ,..., 7 ,^, 
Si,... are given such that 


max 



< min 

l<j<n 



For each k G Z_|_, we set 


— 


1 - 


I3 


^j,k — 


which constitute 


(3.5) 


Pi = {Oj.fc G M : j G {1,.. .,m},k G Z+} , 

P 2 = {cj^k GM:_) G {l,...,n},fcG Z_|_} . 

Note that Pi n P 2 = 0 by (13.5|) . We arrange the elements of Pi and P 2 as follows: 

-Pi = {^0 > 5i > 52 > • ■ ■ } , P 2 = {co < Cl < £2 < • • • } . (3.6) 

For the above parameters, define 

nr=ir(5i+Mn;^ir(i-ci--7i-z) 
nj=n+l r(Ci + 7i^’) ni=m+l r(l - - < 5 ^ 2 ) ■ 
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Note that Pi and P2 are sets of poles of TL^z). To describe the behavior of TLiz) as 
\z\ —>■ 00, we set 


n y m 


i—1 3 — ^ 


and 


+ a;:=j] <5,-^7,, V-= Y[lj t[sf ■ 


i=i i=i 
Due to ( 13 . 3 |) . 

|H(a + i6)r-“-‘''| 

as a —> 00 and 


i=i i=i i=i i=i 


(£)-0\a .,(0,00) (3.7) 


|77(a + i5)r-“-‘^| ^ 

as a —>■ —00. By dSH), it follows that 


“l“l /j.\ “l“l 

■n 


r G ( 0 ,00) ( 3 . 8 ) 


( \ —Loa ^ 

(7) re( 0 ,oo) ( 3 . 9 ) 

as |6| —>• 00. 

The Fox H funciton H()))*(r) (r > 0 ) is defined via Mellin-Barnes type integral in 
the form 


yfi 


(r) := H; 


mn 

Ufl 


:=H 


yfi 


[c, 7 ] 

[ 5 ,( 5 ] 

(ci, 7 i) 
(^i) <5i) 


(Ci/, 7;^) 
(5/i, (5^) 


1 

27ri 


'H{z)r ^dz. 


( 3 . 10 ) 


In (I 3 . 10 p . L is the infinite contour which separates all the poles in Pi to the left 
and all the poles in P2 to the right of L. Precisely, we choose L as follows: 

(i) if Cli > 0 , then L = which is a left loop situated in a horizantal strip 
(or left Hankel contour), runs from —00 + ihi to £ + ihi, and then to l-\-\h2 
and finally terminates at the point —00+i/12 with —oo</ii<0</i2<oo 
and 


5o<f<Co, (3.11) 

(ii) if Ct) < 0, then L = which is a right loop situated in a horizatal strip 
(or right Hankel contour), runs from +oo+i/ii to £+ihi, and then to £+ih2 
and finally terminates at the point +oo+i/i2 with —00 < /ii < 0 < /12 < 00 
and (13.lip . 

(hi) if (V = 0, then L = and L = and for r G (0,7) and r G (i?, 00) 
respectively. 

The following proposition shows that integral p3.10l) is well defined and indepen¬ 
dent of the choice of /ii, /12, and i G (60, Cq)- 
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Proposition 3.1 i |141 Theorem 1.2]). Assume LS. ,5|) and choose the contour L 
as above. Then Mellin-Barnes integral SS. 1 ffj) makes sense and it is an analytic 
function of r G (0,oo) and of r G (O,??) U (? 7 , oo) if uj 0 and uj = 0 respectively. 
Furthermore, 

(i) if w > 0 , then 

OO 

H;7(r) = ^ Res.^s^ [^zy-^] (3.12) 

(ii) if a; < 0 , then 

oo 

= - E [n{z)r-^] (3.13) 

k =0 

(hi) if w = 0, then (|3.12|1 and (|3.13[1 hold for r € (0, 77 ) and r G ( 77 , 00 ) respec¬ 
tively. 


Proof, (i) Let w > 0 and r G (0,oo). Choose hi < 0, ^2 > 0, and £ G R so that 
(13.1111 holds. Take a sufficiently large p G Z+ so that < 0. Then there exists a 
real number M = M{p) > 0 such that 

-ip<M < -ip+i. (3.14) 

Dehne a closed rectangular contour which can be decomposed into four lines 

Cm = Li U L 2 U T 3 U L 4 


where 

Li := {z G C : = i,hi < S[ 2 ;] < / 12 }, 

L 2 := {z G C : K[z] = —M, hi < 3[z] < h 2 } , 

L 3 := {z G C : —M < 5R[z] < £, 3[z] = hi} , 

L 4 := {z G C : —M < 5R[z] < i, 3[z] = / 12 } . 

Note that 'H{z)r~^ is a meromorphic function on z G C \ Pi U P 2 and 



|dz| < 00 . 


Due to (13.8|) (* = 1, 2) 

/ — M p—M 

TL{t + \hi)r~^~'^^*dt < / \TL{t + \hi)i 

-OO j —OO 




dt 


< 




ydt < 00 , 


and 


f 


rh 2 

/ K{z)r~^dz 

= 

/ H{-M + \t)r^-'yt 

Jl2 


Jhi 



rh 2 



\n{-M + it)r^-^*\dt 

hi 





(3.15) 


(3.16) 


(3.17) 
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By (|3.15[l - lj3.17[l . integral (13.1011 absolutely converges and makes sense with L = 

^Ha- 

We next show (j3.12|l . Note that (j3.16|l and (j3.17ll converge to 0 as M —)> oo and 
by the theory of residues, 

^ [ n{z)r-^dz = . 

2^1 Jcm 

Thus our claim follows immediately if we prove the convergence of the residue 
expansion. Let q € N {q > p) he arbitrarily given. Then we can choose a real 
number N = N{q) > 0 satisfying (13.141) where p and M are replaced by q and N 
respectively. Set 

C]\^ := L-^ U L2 U U 

where 

L[ := {z G C : K[z] = -N, hi < 3[z] < / 12 }, 

Lg := {z G C : —N < 5R[z] < —M, = hi} , 

:= {z G C : -A^ < K[z] < -M, 3[z] = /i2} . 


Observe that 


1 

27ri 


f nyy ^dz= Res^^j^ ['H(z)r ^] . 

C T,jr h — 


(3.18) 


k—p-\-l 

By replacing M by iV in (I3.16P and (|3.17|) . 


lim 

N—^OO 


/ n(zy-^dz 

= lim 

[ n{zy-^dz 

Jl{ 

N,M—yoo 



= 0 , 


which implies 


= lim 
N.M—^OO 




r ^dz = 0 . 


lim Y Fies^^j [nizy, 

p,q—¥oo ^ L 

k—p-\-l 

Thus (13.121) is proved. 

(ii) The case w < 0 is an analogue of the case w > 0. By (13.7p . for sufficiently 
large M > 0 (i = 1 , 2 ) 


'H(t + \hiy ^^*dt 


< M 


< 


' M 


|H(t + ihiy 


—t—ihi 


< 


dt 

L © “’( 7 )* 


and 


f n{zy-^dz = l'"n{M + it) 

«/ 2 « /i 1 


.M-ii 


dt 


< 


< 


J \niM + ity^-''*\dt 


(3.20) 


Note that both (13.191) and p.20l) converge to 0 as M —>■ 00 . Then we obtain our 
desired result by replacing hk and M in the proof of the case a; > 0 by Cfe and —M 
respectively. 
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(iii) Finally we consider a; = 0. Note that (13.1611 and (13.1711 hold if r < 77 , 
consequently (j3.12|l follows. If r > 77 , then we have (13.1911 and (j3.20|l which give 
(I3.13|l immediately. The proposition is proved. □ 


In the remainder of this section, we assume 

a* > 0. (3.21) 

Under (j3.21ll . we can choose a contour L = Lsr which is a vertical contour (or 
Bromwich contour) starting at the point £ — ioo and terminating at the point t + ioo 
where i satisfies (I3.11|l . Actually, H()))‘(r) does not depend on the choice of L due 
to the following proposition and it is an analytic function of r S ( 0 , 00 ) (it is 
a holomorphic function of r G C in the sector |argr| < See m Theorem 
1 . 2 .(iii)]). 


Proposition 3.2. Under i3.21\) . Mellin-Barnes integral iS.KA) makes sense with 
L = Lbt- Furthermore, for r G (0, 00 ) 


[ n{z)r-^dz = 7 ^ / Uizy-^dz (3.22) 


ifuj>0, 


1 

27ri 


if uj < 0. If (jj = 0 then 


- f n{z)r-^dz = — f n{z)r-^dz (3.23) 

1 Jlb. 2711 

hold for r < rj and r > rj respectively. 


and 


Proof. We only prove the case w > 0. The proof of the other case is almost same. 
Fix r G (0, 00 ) and 

L = Lsr = { 2 : G C : 5R[2] = f} , 
where I satisfies (13.1111 . Note that the relation 

|H(£ + h)r-^-‘*| - e-‘^?7V-^|t|‘^^+^exp 

holds as |t| —>■ 00 uniformly in I by (13.911 . Hence the contour integral 




absolutely converges and makes sense, 
show that 


n{z)r-^dz 

Due to Proposition 13.11 we only need to 


1 r 

— / n{z)r-^dz = ^ Res.^g^ [n{z)r-^] . (3.24) 

-I 

Given (sufficiently large) p G Z+, we take M as in the proof of Proposition 13.11 
Define a closed rectangular contour which can be decomposed into four lines 

U U Lf 


Lf := {2 G C : ?ft[z] = £, |3[0]| < M} , 

:= {z G C : U[z] = -M, |3[z]| < M} , 
Lf := {2 G C : -M < 5R[z] < £, ^[z] = M} , 
:= {z G C : -M < 5R[z] < £, 3[z] = -M} . 


where 
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Then by the theorem of residues 
1 

27ri 


ICM 




(3.25) 


k=0 


By (13.81) . for any h G 


\TL{t + ih)r 


—t — ih I 


“Id /^\ Id 


\t\J \v 


as t —>■ —oo. Thus 


p pM 

lim / \n{z)r-^\\dz\= lim / \n(-M + dh 


M 

rM 


< lim 
M—¥00 


< lim 2 
M—¥00 


-M 


iii) 


oM 


(s) 


ujM 


M 


M^dt 


= 0 . 




On the other hand, by (|3.9p . 

p p — Af 

lim / \n{z)r-^\\dz\= lim / |'H(t + iM) 

M-^oo Ji^M I ' ' I ' M-^oo Ji' 

r—/ \ t 
< lim / 

M->oo Jf, \r/ 

[ 


dt 


exp < - 


a*MTY 


dt 


= lim exp < — 


-Mir 


M‘ 


riM^ 


re 


dt 


= lim exp < — 


M—¥oo 




/ r]M^ \ _ f riM^ \ 

I re‘^ j \ re^ j 

In ?7 — In r + CLi(ln M — 1) 


since a* > 0, w > 0, and r < rj. Similarly, 

lim f |■H(z)r“^| Idzl = 0. 

M^-oo I ' ^ I ' 

Thus, by taking p —>■ oo in (|3.25p . 

oo 

y R-es.-sy'H(z)r-^l = lim — / n{z)r-^dz 

^ ■' M-s-oo 2711 JrM 

k=0 

. +/ +/ ) 

M^OO 270 Jl-M J^M j 

= lim — f niz)r-^dz = — f n{z)r-^dz. 
M^oo 2711 Ji^M 2711 

Therefore (|3.24|1 holds, and the proposition is proved. 

Proposition 3.3 m (2.2.2)]). 


□ 




lyfi 


[C:7] 
[5, <5] 


_ n 

— I 


[c,7] (0,1) 

(1,1) [5, <5] 
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Proof. Observe that uj and a* of H™(r) and of 


TTm+1 n 


■T* 

[F7] 

(0,1) ' 

/ 

(Ri) 

[D,5] _ 


are same, respectively. Hence is well-defined with the same contour 

used for Therefore the proposition easily follows from (13.11) . (I3.7I) - (I3.9I) . 

and the definition of the Fox H function. □ 


3.2. Algebraic asymptotic expansions of H"™(r) near zero and at infinity. 

Proposition 13.11 gives explicit power and power-logarithmic expansion of H”™(r) 
near zero for w > 0 and at infinity for w < 0. For the reverse case (i.e.) near zero 
for ct) < 0 and at infinity for w > 0, the following asymptotic expansions hold. 

Proposition 3.4. Suppose i3.21\) holds. Then for sufficiently large p G if 

UJ > 0 

H:7(r) = Res,=j, [n{r)r-^] + 0{r-^^), < M < (3.26) 

as r ^ oo, and if oj < Q 

HL7(r) = -Y [n{r)r-^] + 0{r^), -5^ < M < -5^+1 (3.27) 

fe=0 

as r —^ 0. 


Proof. Braaksma [1] proved (13.261) for w > 0. We follow Braaksma’s method to 
prove (|3.27l) for the case w < 0. 

Let O' < 0. Given (sufficiently large) p G take a constant M > 0 satisfy¬ 
ing p.l4l) . Fix h > 0 and let be a right Hankel contour surrounding L'^^. 

Precisely, is a right loop situated in a horizantal strip runs from oo — i/i to 

—M — ih and then to —M + ih and finally terminating at the point oo -I- \h. Let us 
denote by C a closed rectangular contour which encircles Oq ; • ■ • and satisfies 


Then 


f n{z)r-^dz+ f n{z)r-^dz= f n{z)r-^dz. 

Jc J 

= h (-X 

= -^Res^^5^ [n[z)r-^] -h ^ [^niz)r-^dz. 

L. -n J L rr_ 


(3.28) 


Using (13.7p . (13.81) . and (13.211) . and modifying the proof of Proposition 13.21 
1 /■ ^ °° 


k=0 

1 

27ri It-m 


'H(z)r ^dz, 
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where is a Bromwich contour starting at the point —M — ioo and terminating 
at the point —M + ioo. 

We estimate the upper bound of contour integral along ■ Recall (13.91) : 
|■H(-M + i^)r“-‘*| - e““r7-^r^|tr“^+^exp 
as t —>■ oo. Thus, for r < 1, 



1 

27ri 



7i{z)r ^dz 


1 r°° 

\n{-M+ it)r^-^Hdt 

J — OO 


< 





oM 


nM 


I 


\t\ 


—ujM-\-A 


exp 




< r 


M 


The proposition is proved. 


□ 


If w > 0 and n = 0 (i.e.) P 2 = 0, we have the following exponentially asymptotic 
behavior of (see [141 (1.7.13)]). For the proof we refer the reader to [U 

Theorem 4]. 


Proposition 3.5. Assume (13.211) and oj > 0. Then 


lyfi 

as r ^ OO. 


H:7°(r)=0 r(^+^)/‘"exp<^cos 


g* + Ej=m +1 


Ijtjj '* 


4. Asymptotic estimates oe the Fox H function 
Throughout this section we fix 


d e N, as (0,2), /3 e (0,oo). 

For 7 G [0, 00 ), (T G K. and z € <C define 

^ /X_ r(| + 7 + ^z)r(i + z)r(-z) 

r(—7 —/ 3 z)r(l — (7 + az) 


Observe that 

a* = 2 — a, A=^ + 27 + cr 
For each k G Z_|_ we write 

d 

Cpfe := fc, Oi 



+ 7 + A: 


213 — a, T] = a 


02 ,fe := -1 - k. 


Obviously, (13.51) holds, and thus Oi^/c, 02 ,/c and Ci^k constitute Pi and P 2 respectively. 

Remark 4.1. (i) Note that "Ho-, 7 ( 2 ) has removable singularities at z = 0 and at 
z = —lif 7 = 0 and 7 = /3, respectively. Indeed, by (13.11) . 

, _ ^r(f+ /3z)r(i + z)r(i-z) 
r(l —/3z)r(l — (T + gz) 
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Therefore, 'Ha.o{z) has a removable singularity at z = 0. Similarly, 


’H<t,/3(z) = — 


/?r(| + /3 + /3z)r(2 + z)r(-z) 


r(l — P — Pz)T{\ — a + az) 

Thus, 'Ha,i 3 {z) has a removable singularity at z = —1. 

(ii) Assume that /3 G N and 7 = 0, then by (13.21) . 

.o r(f+ /3z)r(i + z) 

na [j{z) = {2 tt) 2 —-/ 

nLir(|-z)r(i-cr + az) 


32+, 


Hence P 2 = 9- 

(iii) If a = 1 and tr = 0, then 

'Ho, 7 ( 2 :) := 


r(| + 7 + /3z)r(-z) 
r(-7-^z) 


Thus, 02 ,fc = 0 for all k € Z+. 
For r G (0, 00 ), we define 

H..7(r) :=H^23 


(1,1) (l-cr,a) 

(f+7j/3) (1,1) (1+7,/3) 




r ^dz. 


Here 


max ( - 1 ,-^ ) < 4 < 1 


(4.1) 


(4.2) 


(4.3) 


L = Lsr = {z G C : K[z] = fo} 
and is chosen to satisfy (13.111) : 

max ^- 1 , ^ 1 < 4 < 0 

if 7 ^ {0,/3}, 

if 7 = 0 , 

max ^- 2 , -1 - < 4 < 0 

if 7 = /3. If a = 1 and cr = 0, then we take Iq such that 

Since (I3.2ip holds (i.e. a < 2), by Propositions 13.11 and 13.21 the value of ]HIo., 7 (t) is 
independent of the choice of £0 as long as it is chosen as above. 

By Proposition 13.41 we obtain the asymptotic behaviors of Ho -,7 (r) at infinity. 


Lemma 4.2. It holds that 


1 . r(f + 7) 

An = -\ 


r(- 7 )r(i-a) 

for r > 1. In particular, if "f G or ct G N, then 

r(f+ 7 + /3) 


nA) = 




r (—7 — /3)r(i — a + a) 


+ 0(r 


r-i +0(r-2) 
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for r > 1. If P G N and 7 = 0 , then there exists a constant c = c{d,a, p,a) > 0 
such that 

= O (exp{-cri/(2/3-«)|^ 


as r —>■ 00 . 


Proof. Observe that 'Hrr.-yiz) has simple poles at z = = Ci_fc = k for all k G Z+. 

By (I3.26p . for sufficiently large p G Z+, 


U 7 (^) = +^Res^^-c, [Ha,^{z)r + 0{r 


1 


fc=2 


= Res^^fc ['Ha,^{z)] r ^ + 0{r 


fe=0 


for r > 1 . Additionally, if 7 G Z_|_ or cr G N, then 


Res^^o [^< 7 . 7 ( 2 )] = lim 


zr(| + 7 + ^z)r(i + z)r(-z) 


z^o y r (—7 — ,5z)r(l — a + az) J 

r(| + 7 + /3^)r(i + z)r(i-z) ^ 
2 ->-o I r(— 7 —/3z)r(l — (T + az) 


= 0 . 


Hence 


T, 7 (r) = Res^=i['Hcr, 7 (z:)]r ^ + 0{r ^) = ( 2 +"! + ^) 


r(— 7 — /3)r(l — cr + a) 




Finally, assume /3 G N and 7 = 0 . By (14.11) . 

( 1 , 1 ) (l-CT,a) 

(f,/3) ( 1 , 1 ) (l,/3) 


1,7,0 (r-) = 

_ (2t)^/ 31/^ 
27ri 


r(f+ ^z)r(i + z) 


L nLi r(| - ^)r(l -a + az) 


-zz dz 


= (27:)^ 


r 

(1 - cr,a) 



(|,/3) (1,1) (l + i,l) •• 

• (1 + ¥> 1 )J 


For the above Fox H function, we have 


7 Q 

a* = 2 — a, A = - + (T + 2/3 —w = 2/3 — a, 77 = a““/3^. 


Note that 


^ < 1, 2/3-a >2-a >0. 

2 2p — a 


Hence by Proposition 13.51 


^1/3+1 


= o ^ 


(A+i)/2/3-a 


exp < cos 


3 + 1 — a 
2/3-a 


77 (2/3 - a) 


-<w) 


l/( 2 / 3 -a)' 


= o 


(exp{-cri/(2/3-«)|) . 


The lemma is proved. 


□ 
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Next we consider the asymptotic behavior of Hcr, 7 (r) at zero. 
Lemma 4.3. It holds that 




d +27 

r~^ 

: 7 < /3 

r\ lnr| 

: 7 = /3 

r 

: 7 > /3 


as r ^ 0. Additionally, if j — [3 £ Z_|_ or ct + a S N, then 

d+2~, , 

- 7 < 2/3 - I 

7 = 2/3-! 






r \ mr\ 
^2 


7 > 2/3 - ! 


as r —^ 0. If a = 1 and cr = 0, then 


d + 27 


as r —^ 0 . 


Proof. Due to (I3.27|) . it is sufficient to compare the order of residues among z = 
t'l.O)3'i,i,3'2 ,o, and h2,i- 

First, let 7 ^ /3 — !■ Then has a simple pole at z = max{c)i_o, 52 ,o}- If 

7 > /3 - then 0i_o < ^ 2,0 so by (13.2711 

1 


as r —>■ 0. Similarly, if 7 < /3 — then 

r n <i + 27 

EIcr,7(r) = Res^^Bi 0 ['Ha,j{z)r~^] + 0{r~^^’°) - 

as r —^ 0 . 

Next, assume 7 = /3 — ! (i.e. Oi^o = 52 ,o)- Then 'Ha,'y{z) has a pole of order 2 at 
2: = 5o = = 52,0 so that 




= lim — ((2 + l)^'H^^( 2 )r 
z-s--! dz ^ ' 


= Resz=_i ['Her,7(2)] + 


I Inrl 


r(|)r(i-u-a)^ 


Then by p.27|) . we obtain the first desired result. 

Now we assume 7 — /3 G Z+ or u + a £ N. Then one can easily see that 
Res 2 =j ,2 g \Ha,-^{z)r~^] = 0. Hence it remains to compare the order of residues at 
-2 = 5 i, 0 j 5 i.i, and 02 ,i- Following the same argument of the above, we obtain the 
additional result. 

Finally, we assume a = 1 and a = 0. Recall (14.21) . and note 02 ,fe = 0 for all 
k £ Z_|_ which implies 


E r 1 d + 27' 

Res^=0i j, [Ho,7(z)r“"'] ~ 

fc =0 

as r —^ 0. The lemma is proved. 


□ 
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For each q € Z+ we define 



Note that HCT° 7 (r) = ]HIcr, 7 (r) and by Proposition 13.31 

(r) = 

holds and (14.51) is well-defined for each q € 1 +. By (13.11) 



r(f+ 7 + / 3 ^)r(i-z)r(i + z) ^^_, 

r (—7 — /3z)r(l — cr -I- az) 


(4.4) 


(4.5) 


and thus 'H^)y{z) does not have a pole at z = 0 if g > 1. Furthermore, has a 
pole at z = — fc — 1 of order at most 2 for each k G Z_|_. Let us denote by 


n ^ k 

Ci,fe '■= k + 1 , 5i fc :=- — -, 02,fe := ~ k 

the new elements of Pi and P 2 for 111^^7(r). Then due to (|3.26|) again, we obtain an 
analogue of Lemma 14.21 


Lemma 4.4. Let q gN. It holds that 

Ht«^(r) ^ r-i 

as r —>■ 00. Additionally, if fl G N and 7 = 0 , then there exists a constant c = 
c(d,a, /3,a,\q\) such that 

H^«^(r) = O (exp 

as r —>■ 00. 


Proof. The proof is similar to the one of Lemma 14.21 The only difference is ik = 
Ci,fe = k + 1 which implies 


Res, 


R(«^(z)r 


(-l)fe-r(f +7 + /3 + /3fc) 


r (—7 — P — Pk)T{l — a + ap ak) 
for each k G Z+. The lemma is proved. 


(fc-bl)«r 


Qr-k-l 


□ 


Lastly, we present another result which is necessary to obtain the upper estimates 
of classical derivatives of p(t, x). Recall a; = 2/3 — a. Let ni, K 2 , ki, and K 2 denote 
constants 

Ki .— ReS^j;— Oj Q [Rct,7(.2^)] ) kv2 .— ReS2— ^2,0 [Rct,7(-2')] 
ki := lim (z - 0i o)^'Hcr, 7 ( 2 ), K 2 := lim ( 2 : - 02 .o)^R<t. 7 ( 2 ) 
which are independent of g > 1. Note that ^2 = 0 if ct -|- a G N. 
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Lemma 4.5. Let q G Z+ 

(i) 


and 7 G [0, c»). 

— 1, 2/3 — then 






as r ^ 0. Additionally, if j = [3 or cr + a G N, then 



d + 27y fo(r““(2’'^^^)) 

2/3 J ^ 0(r^ I Inr|) 


:7r^2/3-|-l 
: 7 = 2^ - i - 1 


as r ^ 0. If a = 1 and a = 0, then 




as r —^ 0. 

If'Y = (3—^, then 



d + 27 


d +27 , d+ 2^+2 

• KiT 2/3 + 0(r 2/3 ) 


d + 27 

~W~ 


9 

• (k 2 Inr + K2 )f + 0(r) 


as r —^ 0. Additionally, // cr + a G N, t/ien 


Hy(r) 



d+ 27 y 

2/3 / ^ |0(r^|lnr|) :/3 = 1 


as r —>■ 0. 

7/7 = ^ — I — 1 , i/ien 

^ R ' (~ ^ 2/3^'^ ) • +<^(^|lnF|) 

as r —>■ 0. Additionally, // tr + a G N, 3/ie?r 

as r —>■ 0 . 

7/7 = 2/3 — then 

as r ^ 0. Additionally, // ct + a G N, 3/ie?r 

^ R ' (~R^) ' + + 

as r —> 0. 


Proof. Due to Proposition 13.IL one can easily see our assertions hold if a; > 0. 
Hence, we assume w < 0. Recall (|3.27l) and take M > 2 satisfying (I3.14|) . 

(i) If 7 ,5 - 7 /3 - I - 1, and 7 y^ 2/3 - then 5i,o y^ O2.0, 5i,i 7^ 02,0, 

and Oi^o 7 ^ 52,i- Hence T/R/z) has simple poles at Oi^o = 02,0 = —1- 
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Therefore, 

2 


— Resz 


+ 0{r^) 


i=i 


= - {nii^i.oYr + K2{^2fiYr ®^’“) 


- Resz=_2 


n\?hz)r-^] +0(r“”(2’'^^)) 


d + 27 V d +27 


2/3 


/Cir 23 — (—1)'^ • K2r 


(4.6) 


- Res2=_2 




as r —?> 0. Note that ^2 = 0 if 7 = /3 or cr + a G N so the second term of (14.61) 
vanishes. In that case, observe that 


Res„-2|«“ (^)r-l = 


:T# 23-|-1 

lnr|) :7 = 2/3-|-1. 


Furthermore, if a = 1 and u = 0, then 

Res2=r,, = 0 

for all k G Z+. Thus (i) is proved. 

(ii) Suppose j = f3 — Then ho = 5i,o = h 2 , 0 j = K 2 , and ki = /t 2 . Note that 

7i^a)y{z) has a pole at 

» d + 27 
5o =-— = —1 


2^ 

d 


of order 2. By (14.4p . 

Res2=_i = lim ^{{z + lfn,^^iz)z<^r-^) 

L ’' J 2->--i dz ' ^ 

= (-1)'^ •/t2rInr + Res2=_i 'Ui^lyiz) 

= (- 
= (- 

Therefore, 


^ • k 2 r In r + ( 

-1)^ 

• (k2 - qk2) 

r 

' • (k 2 Inr + K2)r - 

- g(-l)‘^d2r 

(4.' 

'-iyk2r + 1 

O(^-in(2,±t0±2), 

1 :/3^1 

0(r 

|lnr|) 

:/3=l 

a G N, then 

k2 = 

0. Thus (ii) 

is proved. 

5i,i = 52,0 a-nd TL 

^)yiz) has a 

pole at 

5i,i — 52,0 - 

■- -1 




(hi) Now we let 7 = /3 — I — 1. The 

d 27 “t“ 2 

W 

of order 2 , and (j4.7l) holds. Also note that H^Jy^z) has a simple pole at 

d + 27 1 


^ — hi^o — — 




2/3 /3 
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Therefore, 

= -Res._ d +27 ni^}y{z)r 

d+2^\'^ d+2~, 


2/3 


- Res2=_i 
Kir 2 ^ + 0 (r|lnr|) 


+ o{r) 


as r —>■ 0 . If (T + a G N, then K 2 = 0 and 


H“ M = - (-^)‘ 


d+2-y 

Kir 2 /a + 0{r) 


Thus (iii) is proved. 


(iv) Finally, we assume 7 = 2/3 — Then Oi^o = 52,1 and Ti^aJiiz) has a pole at 


d + 27 


= 5i.o = 52,1 = —2 


of order 2. By (14.41) . 

Res,=_2 






(-^^) • kir^ Inr + ^ ((z + 2YH.^,{z)z^) r^ 


d + 27 \^ //. 1 \ 2 ( d + 27 

' • (/Cl Inr + Ki) + g - 


2/3 


2P 


q-l 


kir^. 


Note that 'H^}y{z) has a simple pole at z = 02,0 = ~1- Therefore, 


M.^J‘}y{r) = -Res._ d+ 2 ^ 'Hi‘‘}y{z)r 
d + 27 ^« 


- Resz=_i 
,2 


ni^hzy-^ +o(r 2 ) 


2/3 

d + 27 ^'^ 
2/3 


• (di Inr + /Cl) — /C 2 r + O(r^) 

• (/Cl In r + /Cl) + 0 (r) 


(4.8) 


as r —>■ 0. Additionally, if cr + a G N, then /C 2 = 0 in (14.81) hence we obtain the 
desired result. The lemma is proved. □ 


5. Explicit representation of p{t, x) and its derivatives 
T ake the function ]HlCT, 7 (r) from (14.311 . and define 

Write 

pit, x) := po,oit, x) = ■ (5-1) 

Note that ]HIcr. 7 (r) is a bounded function of r G (0, 00 ). As a corollary of Lemmas 
14.21 and [4. 31 we obtain the following upper estimates of Pa-,-yit, x) for 7 G [0, 00 ) and 
(T G K.. 
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Theorem 5.1. Let a G (0,2), /3 G (0,oo), 7 G [0,oo), and ct G M. 
|^|2/3^-a > 2 

\Pa,j{t,x)\ < 

and for \x\'^^t~‘^ < 1 


\P<7,j{t,x)\ < 


r^- 


a{d+2.~,) 


|^|_d_2T,+2/3^ 

|^|_d_2T,+2/3^ 


— <7—0; 

-<7-0 


(1 + \\n{\x\'^H “) I) 


7 < /3 
7 = /3 
7 > /3 


Furthermore, 

(i) If a + a € N, then for \x\'^^t~°‘ < 1 

\Pa,'f{t,x)\ < < (1 + |ln I) 

[|a;|-‘^-27+4/3i-a-2a 

(ii) If a = 1 and a = 0, then for \x\‘^^t~°‘ < 1 


7 <2/3 
7 = 2/3 
7 > 2/3 


Then for 


d 

2 

d 

2 

d 

2 ■ 


d 

2 

d 

2 

d 

2 


d+2-Y 

bo.7(^.a^)l ^ ^ • 

(^m,) If j £ Z+, t/ien /or > 1 

(iv) If fl £ N and 7 = 0 , then there exists a constant c = c{a,P, tr) > 0 such that 

for > 1 , 

\Pap{t,x)\ < \x\-H-^ e^p[-c{t-‘^\x\^^)^] . 

(v) If ^ = 13, then for > 1 

Ux\-‘^-'^l^t-'^+°‘ :ctGN 


\Pcr,p{t,x)\ < 


|^|_d-2/3^-<T :crGK\N 


and for “ < 1 


» —(7 —O— ■: 


: - < /3 

\PaAt^^)\ ^ \ (1 + |lii(|a;p/5i-“)|) : I = /3 

,^,-d+2p^-a-2a : I > /3. 


Remark 5.2. (i) By Theorem 15.11 Pa^^{t,-) £ Li(K‘^) for all 7 G [0,/3] and cr G K. 
Furthermore, if a = 1 and u = 0, then po,~i{t, •) G for all 7 G [0,oo). 

(ii) Observe that 


_ a(d + 27 ) _ Q 

p„A^^x)=t ^x) 


(5.2) 


which implies 

rT 


/ / |pCT,7(t,a::)! 

do aRd 




/ Ipo- 7(1, a:)| da; I dt < c» 
dRa ’ / 


if O’ + ^ < 1. Thus, under this condition, one can consider Riemann-Liouville 
fractional integral and the Fourier-Laplace transform of Pa^^{t,x). However we do 
not use such transforms in this article. 


The following theorem handles the interchangeability of A'^' and D/. 
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Theorem 5.3. Let 7 S [0,c»). For any a G M., m G N, and rj G {— 00 , 1), 

—P(7+m,7 (^7 ; 

and 

-P? 7,7 Pcr+7],^{i: • 

Proof. First we show 


dt 


for any cr G R. By (13.11) . 


az 


—Pa,^it,x) =Pa+l.jit,x) 


1 (7 

+ 


(5.3) 


r(l — cr + az) r(—(T + az) r(l — (T + az) 

which implies 

LLcr^'yi^Z^GZ — ^CT_|_x^7(z) ~t" (j'htfj^.yi^zf 

Then by Proposition 13.31 


2/3 


dt^ V 22/5 


t 


t “ ) = — / 'Ha,y{z)az ( “ ) dz 


27ri 

27ri 




• {'H.i+a,^{z) + aFLa.^i^z)) 


= t 


^cr+ 1,7 




22/3 


.12/3 


Therefore, 






227 


rj^djl 

= Pa+\aif,x) 


d- 27 ^-a-l]] 


l-(7+l,7 


22/3 

W 


‘■tT,7 




22/3 


t 


and (15.31) is proved. Thus to complete the proof of the theorem, due to (12.21) and 
(lOl) . it is sufficient to show that for tr < 0 and 77 < 1 

~ Pri+<T,^itT x). 

Take io > —satisfying (I3.11|) . For 5i[z] > —it holds that 

[\t - s)-^-h-^+^^ds = r(l-i7 + az) 

Jo T{l-a-p + az) 


n-a) 




r(l — 77 + az) 
r(l — (T — 77 + az) 


— ^77+cr,7 ( 2 ^) • 


Observe that 
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By (13.91) and the Fubini theorem, 

= /" . — p^^y{s,x)ds 

Jo i 

2^^TT-i\x\-‘^-^^ /■* r (t-s) 


27ri 

2 ^'>Tr-i\x\-‘^-‘^'^ 

27ri 

2 ^'^TT-i\x\-‘^-‘^'^ 


0 JL 


r(-(T) " [ 22/3 


n|2/3< 


dzds 


n fz) (^ 


/' 

VO 


(t- s)- 


27ri 






22/3 


r(-cT) 


dz 


-s-’/+“^ds 


dz 


= 22'/7r-2|a;r‘^-2^t-'' 
— Pr)+cr,-i{t-, x). 

The theorem is proved. 


l^?7+cr,7 




22/3 


□ 


Remark 5.4. Let 0 < e < T. Then by Theorems 15.11 and 15.31 p{t,x) and ^{t,x) 
are integrable in cc G uniformly in t G \£,T]- Thus 

e~^^^’^^p(t,x)dx= I e~'^^’^'>^{t,x)dx. 
at j^d j^d at 

To estimate the classical derivatives of p{t,x), we need the following theorem. 

Theorem 5.5. Let a,/3, 7,17 be given as in Theorem 15. f I and n G N. Then for 

\x^2Pt-a > I 


and for “ < 1 


exp |-(|a:|2/3t-“) 1 

: /3 G N ,7 = 0 

- 27 —2/3—n^—i7+a 

: 7 G or (T G N 

- 27 —n^—cr 

: otherwise, 

^ _ a(d + 2 ~^+ 2 ) 


2-n-f. <y - 2 ^ - 

+ |ln|x|2/3t-“|) 

:7 = /3-i-l 

— (i— 27+2/3—n^—fT—a 

:7>/3-i-l. 


Additionally, for < 1, the followings hold: 

(i) If a = I and ct = 0, 


\D:poAt,x)\ < 


(ii) If a + a gN, 


2-nj.-o-- 




»(d+2^ + 2) 

W 


|£i>^,.y(Lx)| < <( |a;|2-'^t-<^-2«(i + |ln|a;|2/3t-“|) 
1 ^ 1 —d— 27 + 4/3 —n^ — cr—2a 


7 < 2/3 - I - 1 
7 = 2 ^ - f - 1 
7 > 2 ^- 1 - 1 . 
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(in) If-f = p, 


|^|2-n^-T-a-^ . I + 1 < ^ 

\DOp,^ 0 it,x)\ < <( |a;|2-"t—2“(1+|ln|a;|2/3t-“|) : | + 1 =/3 

. I + 1 > 

Proof. Write R = R{t,x) := 2~‘^^\x\‘^^t~°' and recall from (14.511 . By Proposi¬ 
tion [2i31 and (13.11) . 

(i?) = -2/3 ^hS^+i) (i?) 


y{d+2) 


for each q G and z = 1,..., d. Hence 

\DooiPan{t,x)\ 

227 


rd /2 


^i|2.|-d-27-2^-<. ^ 27)H,,^ (i?) + 2/3HW (R) 


< C\2 

and 

227 


I—d —27 —Ij.— cr 


{d + 27 )H,,^ (i?) + 2/3H« (i?) 


027 


rd /2 


(^(^ + 27)H.,7 (i?) + 2/3HW (i?) 

- 5,j\x\-'^-^^-^ {(d + 27)H,,^(i?) + 2 ^hW (i?)| 

+ (d + 27 + 2)x^x^\x\-‘^-^^-^ {(d + 27)H,,^(i?) + 2/3HW (i?)} 
+ 2Px^x^\x\-<^-^^-^ {(d + 27)HW (i?) + 

2 

< |(d + 27)H^«-i)(i?) + 2/3H^«^(i?) . 

q=l 

Inductively, for any n £ N, 

n 

\DZp.,^{t,x)\ < |(d-f 27)Hi';-i)(i?) + 2/IHi«) (i?)| . (5.4) 

9=1 

By Lemma 14.2114.41 and (15.4p , 

f exp {-c(|x|2dt-“) 23^ J 

\Dlp<y,^(t,x)\ < \ | 2 .|-d- 27 - 2 / 3 -n^-a-Ha 
[^| 2 .|-d- 27 -n^-<T 

for > 1. To estimate the upper bounds for \x\^^t~°‘ < 1, observe that 


7 = 0 , /3 £ N 
7 £ Z+ or (T £ N 
otherwise 


(d -f 27 ) - 


d -I- 27 

w~ 


9-1 


+ 2/3 - 


d -I- 27 


= 0 . 





















ASYMPTOTIC BEHAVIORS OF FUNDAMENTAL SOLUTIONS 


27 


Thus, by Lemma [4.51 and (15.411 . 


(d+27 + 2) 




| 2 -lrtU-o --2 

\D'^Pa,^{t,x)\ < { “(1 +|ln|a;|2/5r“|) 

|^|—d—27+2/3—|n|^—IT—a 

for R<1. 

Additionally, if cr + a £ N, by Lemma 14.51 

\Dlp^^^{t,x)\ < { |a;|2-l"lf-'^-“(l+|ln|a;|2/5i-«|) 

|2.|-d-27+4/3-|n|^-cr-2a 


for i? < 1. Also, if 7 = /3, 


. | 2 -„^-a-a-. 


»(rf+2) 

2/3 


7<^-|-l 


7 < 2/3 - f - 1 
7 = 2/3 - I - 1 
7 > 2/3 - I - 1 


: f+ 1< 


\D^Po^p{t,x)\ < \ |a;|2-l"lf-'"-2“(l +|ln|a:|2/5t-“|) ; d + l = ^ 

: f+ 1>^ 


|2.|-d+2/3-|n|^-CT-2a 

for i? < 1. The theorem is proved. 


□ 


6 . Proofs of Theorems 12 .1 L 1^75] and I^TTl 


By Remark !?/^ Pa,-yit, ') G if 7 G [0, /3] and cr G M, andpo, 7 (i; •) € Li(K‘*) 

if a = 1, 7 G [0,oo), and cr = 0. Due to Lemmas 14.2114.31 Theorems 15.11 1??^ and 
15.51 to prove our desired results, it is enough to show 

-^{p..7(iG)} = (e.i) 

which is equivalent to 


Pa,^{t,x) = A'^p^^o{t,x) = A'^I}^p{t,x). 
We devide the proof into the cases d> 2 and d = 1 . 


Case 1 : d> 2. We choose £0 G R in (14.3p such that 


if 7 = /3, 


if 7 G (0,/3), 


max(-2, -1 - < 4 < -1 




/ 1 d-1, 

max(-l,-,-— < 4 < 0 

a Afj 


if ^ = 0. If a = 1, 7 7 ^ 0, and <7 = 0, then we take io such that 

^ p ^ -1 

/3 4/3 ° p- 


( 6 . 2 ) 


(6.3) 


(6.4) 


(6.5) 


Under the above restriction on (Ii31) is well-defined and the value of Po -,7 (tj x) is 
independent of the choice of £o- 
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Due to the Fourier transform for radial function (see |341 Theorem IV.3.3]) 




2 #+27 


lei^ 




where Jd_-^ is the Bessel function of the first kind of order | — 1, (i.e.) for r € [0, oo), 

(- 1 )'= 






k!r(k + l) 


(d 


2fc-l+f 


It is well-known if to > —1 then 

Jm(t) = 

Due to (j6.2p - (j6.5l) . 


Oit"') :t^0+ 
: t —>■ oo. 


^OO ^1 /'OO 

/ dp< p-i-^^-^0^°-^dp<oo. 

Jo ^ Jo Jl 

Recall 5R[2:] = £o along L. By the above and (13.9F 

p2P^-a' 


10 JL 


Therefore, by the Fubini theorem 


22/3 




\dz\dp < oo. 




27ri 




J4_i(|^|p)dp 






By using the formula [TJ (11.4.16)], 

r p-i-27-2/3. i(|C|p)dp = 2-i-27-2/3.|^|#+27+2/3.-l 

Jo 1^1 r(| + 7 + ^z 


/X r(- 7 -^z) _ r(z-h i)r(-z) 
«-. 7 l^lr(f +7 + / 3 z) r(l-a + az)' 


we have 
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Hence 


27ri 




/i-«\ ' 


f r(^ + i)r(-z) 




27ri 


II r(l -cr + az) 


^^, 7 ( 2 ) ^ 22/3 j 






r(i-^)r(z) 
i r(l - cr - az) 


(|e| 2 ^i“) "dz 


( 0 , 1 ) 

(0,1) (cr,a) 


Therefore, 


-^{p..7(ic)} = |ei"'^i-"i?a,i-.(-|?P''i“)- 


Case 2: d = 1,7 G (0,/3). We choose £0 such that 

7 1 

max(-l,-- - —) < 4 < 0. 
Since Pa.-yit, x) is an even function, 

/OO 

e~"^'"pa,~f{t,x)dx 

-00 

poo 

= 2 / p^ .y{t,x) cos{^x)dx 

/o 


0 


2 / Per 7 (t, a:) cos(^a;)dx + / po- 7 (t, x) cos(^x)da: 

\ Jo ’ Jt“/2/3 ’ y 

‘1 ^00 

PCT, 7 (l,a:) cos(t^^a;)dx + / p„(1, x) cos{E^^x)d:i 


= 2r'^—f 

The last equality holds due to (15.21) . Set Bromwich contours 

L<c := {z G C : 5J[z] = £<c}, T>c := {z G C : 3?[z] = £>c} 

where 

max(-l,-^ - ^) < <-^, -^<J>c<0. 


By (I3H), 


0 Jl^ 




- 27 - 2 / 32-11 


|dz|da 


+ f f \2^^^'Ha^^{z)x ^1 |(iz|da; < 00 . 

Jl JL>e 


Therefore, by the Fubini theorem. 
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/ Pcr,'y{^,x)cOs{t'^ ^x)dx 

Jo 


227^-1/2 rl 

'0 JL^c 


27ri 

227,1-1/2 


2‘^<^^n^^^iz)x-^'<-^f^^-^cos{t^Cx)dzdx 


27ri 


[ 2^^^Ha,j{z) [ cos(t^S,x)dx 

L/ cr r 0 


dz. 


Similarly, 


/ CO 

Pa,'y(^, x) COs{t^ ^x)dx 

2277 ,- 1/2 r „ r 

■ / 2^^^H,,^{z) / 

J L-^c .*/ 1 


27 2/3z ^ nnaf'tJe 


27ri 


cos(t^^ ^x)dx 


dz. 


For 77 G R, it holds that (see [301 (1.8.1.1), (1.8.1.2)]) 


/ 


X ^ cos{r]x)dx =-- 


2X 


, 5?[A] < 0 


and 


/ 


“ 2 A 1 . r(-2A)cosA7r ii ^2 (-A; i, 1 - A;j 

cos{r]x)dx = — - - - ' -^^ 


5i[A] >--, 


|r;|-2^ 2A 

where 1 F 2 (~A; i, 1 — A; denotes the general hypergeometric function, i.e. 

r(p + fc) 


^ , r(g)r(r) ^ 

iF 2 (p; q, r; z) = 2_^ 


^g^id + k)T{r + k) k\ 


, pG C, g,r G C\{0,-l,-2,---}. 


Observe that for 7 + /3z G C \ N, 


ii"2 (-7-/9^;^,1-7-/3^;-^^ 


0rr(l - 7 - ,8z) r(-7-/3z + fc) ( /^) 

tTT 


n-j-^z) ^^T{^ + k)r{l-^-/3z + k) 


k\ 


= -E 


y7r(7 + /lz) 




fc^or(| + fc)fc!(-7-/3z + fc) 


= i-E 


^ + (dz 


_ 

^ (2fc — l)!2fc (—7 —/3z + A:) y 4 


=: Ft,^iz)- 


Since the series in Ft^^{z) converges absolutely for fixed t, one can easily see that 
Ft^^{z) is holomorphic in {z G C : < K[z] < —^ 2 ^} by Morera’s theorem. 

So we obtain 

'■i ^ 22^71-1/2 

Pct. 7 ( 1 , x) COs(t 2/3 ^x)dx =-—^ 

2711 
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and 


/ OO 

Pcr,7{lj x) COs{t^ ^x)dx 
027 - 1/2 / f 

= —- / 22 '^'"H,T, 7 ( 2 :)r(- 27 - 2 / 3 z)cos ((7 + /3z)7r) 

\ ^>c 




27 + 2/32 


dz 


If 7 ^ /3, by (123]), 


T~Lcr,^{z) r(i + 7 + ^z)r(i + z)r(—z) 


2(7 + /3z) 2 r(l — 7 — /3z)r(l — cr + az) 

Thus z = — ■? is a removable singularity. These facts lead to 


Res,^_: 


22/32 R+7(^) / n 

, 2(7 + /3z) 


= 0 


( 6 . 6 ) 


and 


27ri 

Therefore, 


L /■ Ft^{z)dz. 

TTI 2(7 + / 3 z) ’ 270 Jl 2(7 + ^z) ’ 


poo 

/ Pct. 7(1, ic) cos(t^^a;)dx 

io 

9277,- —1/22^ I c|27 r 

= -- / 2 ^^^'Hcr. 7 (z)r(- 27 - 2 / 3 z)cos ((7 + /3z)7r) (t“|^P^)^dz. 

L'-tr- 


By (1331), 


'H<,, 7 (z)r(- 27 -2/3z)cos((7 + /3z)7r) = 71^22 

1 (1 — cr + az) 


Hence 




poo 

{Pa, 7 (i,-)} (6 = 2 t"'^"^ / p„, 7 (l,a;)cos(t^^a;)(ia; 

Jo 


27ri 


, r(l - cr + az) 




Case 3: d = 1 ,7 = 0. Let 

max(-l,--,--;^) < 4 < 1 . 
a zp 

Again, we follow the argument in Case 2. Note that 

n.,o{z) r(i + /3z)r(i + z)r(-z) 

2/3z 2 r(l - /3z)r(l - cr + az) 

has a removable singularity at z = 0 if and only if cr G N. Thus (| 6 . 6 I) holds if cr = 1 
and we immediately obtain 

-F{pi.o(i, •)} = t-^E^,o{-\e^n = 



















ASYMPTOTIC BEHAVIORS OF FUNDAMENTAL SOLUTIONS 


32 


By Remark 15.41 and Theorem 15.31 

•)} = •)} = •)} ■ 

Thus, 

T{pit,-)}{o = Eu-\e^n+Rio- 

By dOl), 

E{p{t,-)} {^) = f e~'^^p{t,x)dx = f x)dx 

Jr Jr 

Jr 

Then by the Riemann-Lebesgue lemma, -7^ {p(t, •)} (C) converges to 0 as t —>■ oo. 
This implies R{^) = 0 and we obtain the desired result. 


Case 4: d = 1, 7 = /3. Now we additionally assume cr + a £ N and take £0 so that 

max(- 2 , “1 “ ^) < -^0 < 0 . 

Note that every argument in Case 2 holds except that 

n.A^) r(i + /3 + /3z)r(2 + z)r(-z) 

2/3(1 + z) 2r(l — /3 — /3z)r(l — a + az){z + 1) 
has a removable singularity at 0 = — 1 if ct + a £ N. Thus (| 6 . 6 p holds and we obtain 
(EH). The theorems are proved. 
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